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2 ■ ABSTRACT 

' A number of important processes taking place around strong shocks in supernova 

d . remnants (SNRs) depend on the shock obhquity. The measured synchrotron flux is a 

function of the aspect angle between interstellar magnetic field (ISMF) and the line 
of sight. Thus a model of non-thermal emission from SNRs should account for the 
^ . orientation of the ambient magnetic field. We develop a new method for the estima- 

' tion of the aspect angle, based on the comparison between observed and synthesized 

radio maps of SNRs, making different assumptions about the dependence of electron 
fT^ I injection efficiency on the shock obliquity. The method uses the azimuthal profile of 

CN| ' radio surface brightness as a probe for orientation of ambient magnetic field because 

it is almost insensitive to the downstream distribution of magnetic field and emitting 
electrons. We apply our method to a new radio image of SN 1006 produced on the 
' basis of archival VLA and Parkes data. The image recovers emission from all spatial 

. structures with angular scales from few arcseconds to 15 arcmin. We explore different 

models of injection efficiency and find the following best-fitting values for the aspect 
. ^ , angle of SN 1006: 0o = 70° ±4.2° if the injection is isotropic, (/-o = 64° ± 2.8° for quasi- 

. perpendicular injection (SNR has an equatorial belt in both cases) and 4>o = 11°±0.8° 

' for quasi-parallel injection (polar-cap model of SNR). In the last case, SN 1006 is ex- 

\ pected to have a centrally-peaked morphology contrary to what is observed. Therefore, 

our analysis provides some indication against the quasi-parallel injection model. 

Key words: ISM: supernova remnants - shock waves - ISM: cosmic rays - radiation 
mechanisms: non-thermal - acceleration of particles - ISM: individual objects: SN 1006 



1 INTRODUCTION 

Non-thermal emission of supernova remnants (SNRs) is in- 
tensively studied. It carries out important information about 
physics of strong shocks, kinetics of cosmic rays and mag- 
netic field properties. SNRs are therefore observed with 
space and ground based observatories in the X-ray, 7-ray 
and radio bands, and modeled with advanced codes for MHD 
and /or particle kinetics. 

The interstellar magnetic field (ISMF) creates different 
obliquity angles with the shock normal in different places 
of the SNR surface. Efficiencies of injection and accelera- 
tion, compression and/or amplification of ISMF may de- 
pend on the shock obliquity. Therefore, in order to model 



non-thermal emission, it is important to make assumptions 
about ISMF orientation around these objects. 

In many cases the ISMF, galactic in origin, may be 
assumed to be rather uniform on the scales of SNR sizes. 
This is likely to be the cas e of the bilateral su p ernov a 
remnants (here after BS NRs, Kesteven fc Caswelll l|l987h : 
iGaensled (119981 ): but see lOrlando et all |20o3) for impUca- 
tions of asymmetries in BSNRs), which are characterized by 
strong, opposite limbs. BSNRs with symmetric structure in 
radio images limit the orientation of the ISMF component 
in the plane of the sky. In the case of the BSNR archetype 
SN 1006, for instance, the ISMF may be parallel to the sym- 
metry axis, spanning from SE to NW, or be perpendicular 
to it, running from NE to SW. The former corresponds to 
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barrel-like (eq uatorial belt) structure and the latter to polar- 
cap structure (|Rothenflug et al.ll2004l V However, in order to 
draw reliable conclusions about the shock obliquity, it is also 
necessary to consider the aspect angle - the angle between 
ISMF and the line of sight - which is still unknown. 

In the present paper, we propose a new method for de- 
termination of the aspect angle. We extract the azimuthal 
profiles of the synchrotron surface brightness distribution in 
a given SNR and compare the observed profiles with those 
synthesized from theoretical models, making different as- 
sumptions on the aspect angle and/or on the details of in- 
jection and acceleration of electrons. The "true" aspect angle 
is that of the best-fitting model. As a first application of the 
method, we analyzed SN 1006. 

The plan of the paper is as follows. The method itself is 
explained in Sect. [2] where results of numerical simulations 
are presented. In Sect.[3]we use our method to determine the 
aspect angle toward SN 1006. To achieve this, we produced 
a new radio image of the remnant, obtained from re-analysis 
of archival VLA and Parkes data. An approximate formula 
for the azimuthal variation of the radio surface brightness is 
presented in Sect. |3] where we also discuss effects of the shock 
modification on our results as well as some consequences on 
the injection model in SN 1006. 



2 ASPECT ANGLE IN BILATERAL SNRS: 
THE METHOD 

The method to determine the aspect angle (jjo between ISMF 
and the line of sight on the basis of radio emission in SNRs, 
can be deduced from the simulated radio maps produced for 
a Sedov SNR expanding through a uniform ISM and ISMF, 
presented in Fig. [1] As shown in the figure, the azimuthal 
profile of the surface brightness is sensitive to the aspect 
angle: it is constant for (po ~ 0° (Fig. [T^: ISMF is directed 
towards the observer) and it is steepest for (po ~ 90° (Fig.[lJ;: 
ISMF is in the plane of the sky). Thus, comparison of an 
observed azimuthal profile with theoretical ones allows one 
to conclude about the aspect angle. 

The radio surface brightness S at some 'point' of SNR 
image is 

S<x J KB^'+^^^^'dl (1) 

where s and K are the index and the normalization of the 
electron energy spectrum, respectively, B is the magnetic 
field (MF) strength and the integration is along the line of 
sight within the volume of SNR. For a Sedov SNR in a uni- 
form ISM and uniform ISMF, the downstream distributions 
of K and B are self-similar, i.e. may be written in the form 
K = K^{Qo)K{f), B = Ba(eo)B(f , Qo) where r is the radial 
distance from the centre of SNR, Oo is the shock obliquity 
angle between ISMF and the shock normal, index 's' marks 
the immediate post-shock values and the upper bar marks 
the downstream variables normalized at their own values at 
the shock front (|Revnoldslll9"9i : |Petrukll2006l ) . 

The post-shock magnetic field is generally a subject of 
compression, and under conditions of efficient cosmic ray 
acceleration, amplification. The magnetic field can be ex- 
pressed as Bs = ylB(So)-Bo where Ab is a product of the 
compression factor (jb and an amplification factor. Since it 



is unknown whether magnetic field amplification depends 
on the shock obliquity, we assume it to be independent of 
Oo. In this case, the variation of the post-shock MF with 
obliquity is only determined by the compression l|Revnoldsl 
\l9m . Bs cx aB(eo)Bo, where 

/l^o^tan^y^' 
^ 1+tan^eo ) ' 

a — 4 is the shock compression ratio for unmodified shocks 
(changes in this prescription and our results due to the shock 
modification are discussed in Sect.|3)). 

At the shock, the normalization, oc ^(Oo) where ? is 
the injection efficiency defined as the fraction of accelerated 
electrons. There are three alternatives for dependence of in- 
jection efficiency ? on obliquity of the shock typically consid- 
ered in the literature: isotropic injection (i.e. <; independent 
of &o), quasi-parallel ( ■; oc cos^ 8s) or quasi-perpe ndicular 
(? oc sin^ Gs) injection (jFulbright fc Revnoldj|l990l ). There- 
fore, the injection efficiency either decreases (quasi-parallel) 
or increases (quasi-perpendicular) with increasing obliquity, 
or it is independent of Oo. The MF compression factor in- 
creases with Oo. 

In symmetric, bilateral SNRs (like SN 100©, the pos- 
sible orientations of the ISMF in the plane of the sky are 
limited. Namely, ISMF may be parallel to the symmetry 
axis or perpendicular to it. For example, in SN 1006, if the 
injection is isotropic or quasi-perpendicular, the bright limbs 
correspond to the magnetic "equator" (equatorial belt) and 
ISMF should be oriented in SE-NW direction. If injection 
prefers quasi-parallel shocks, the bright limbs of SN 1006 are 
two polar caps and MF should be oriented in the NE-SW 
direction. In other words, the model of injection determines 
the orientation of the plane-of-the-sky component of ISMF. 

There are theoretical expectations that i njection is 
higher at parallel shocks l|Ellison et all 1 19951 : IVolk et al.l 
[2003). In contrast, observational evidence seems to ar- 
gue for isotropic or quasi-perpendicular inj ection in SNRs 
(jFulbright fc Reynold j|l"990l : lOrlando et al.ii2007) . We con- 
sider all three possibilities in this work. 

We synthesized a number of radio surface brightness 
maps of a Sedov SNR evolving i n a uniform ISM and uniform 
ISMF. We used the m odel o f Rcynoldj l|l998l ) with some 
extentions presented in IPetruk (2006,) . This model is able 
to account for the evolution of Ks oc V~'' and Ab oc V^/^ 
where b and l3 are constant and V is the shock velocity. 
These relations reflect an eventual evolution of the injection 
efficiency and MF amplification at the shock. 

Radio maps of the model SNR were calculated for a 
range of indices: s from 2 to 2.2, b from —3/2 to 2, P from 
to 2. For each set of parameters, we produced a series of 
images by changing the aspect angle tpo from to n/2. For 
the image with (^o = 7i"/2 in each series, we found the radius 
of projection £>smax which corresponds to the position of the 
maximum brightness Smax = 5(£>smax, ^Smax) (by definition 
</5Smax — Tv/2 for tMs positioujfl- We then traced the surface 

^ Note that SN 1006 is symmetric with respect to the axis be- 
tween the two radio lobes. However, the lo b es app ear slanted and 
converging to the SE (see lOrlando et al.l ||2007| ). for a possible 
explanation of this kind of asymmetry). 

^ The radius f)Smax is almost the same for any aspect angle, ex- 
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Figure 1. Surface brightness distributions (in linear scale) of the radio emission of adiabatic SNR for different aspect angles: if'o = 0° 
(a), 60° (b), 90° (c). Uniform ISM, uniform ISMF and isotropic injection are assumed. In cases (b) and (c), the component of ISMF in 
the plane of the sky is parallel to the horizontal axis. In all three plots, the levels of brightness are spaced in the same way. 



brightness S(£ismax,V') at this fpsmax for azimuthal angles tp 
from to 7r/2 and plotted these distributions normalized by 
the brightness S(£ismax, i/JSmax)- 

Fig.[2]shows the plots for different aspect angles, assum- 
ing & = 0, /3 = 0, s = 2.2, and considering the three models 
for obliquity dependence of injection efficiency. We found a 
reasonable result, namely that all the azimuthal profiles of 
surface brightness which we obtained are almost insensitive 
(within 10 — 20%) to the values of s, b and (3 (at least in the 
case of uniform ISM and uniform ISMF assumed here) . That 
means the azimuthal profiles of radio brightness are almost 
independent of the shapes of the distributions of relativistic 
electrons and magnetic field downstream of the shock. Such 
stability allows one to safely use the method proposed here 
for determination of the aspect angles from radio maps of 
SNRs if the ISM and ISMF in which SNR expands can be 
considered to be mostly uniform. 



3 THE ASPECT ANGLE IN SN 1006 

The method is applicable to SNRs expanding in almost uni- 
form ISM and ISMF. BSNRs, particularly SN 1006, are ideal 
candidate targets of this study. In this SNR, uniform condi- 
tions are likely to be achievable since SN 1006 is located over 
500 pc above the galactic plane. In addition, the SE edge of 
the SNR exhibits a near-spherical shape, a good argument 
for expansion of the shock into a uniform ISM. 

3.1 Radio data 

We now apply our method for the determination of the as- 
pect angle to SN 1006 considerly the best example of a sym- 
metric, bilateral SNR. To this end we produced a new radio 
image of SN 1006 at A ~ 20 cm. 



This image was produced on the basis of archival VLA[f| 
data obtained in October 1991, February 1992 and July 
1992 in the hybrid AnB, BnC and CnD configurations, re- 
spectively. The observations in the AnB configuration were 
carried out at 1370 and 1376 MHz, while the observations 
in the BnC and CnD arrays were performed at 1370 and 
1665 MHz. The data corresponding to the more compact 
configurations of the VLA, BnC and CnD, were published 
as a part of an expansion study of SN 1006 (iMoffett et al.l 
Il993h . but not the data from the AnB configuration, which 
provides the highest angular resolution of southern sources. 
The new interferometric image is produced on the basis of 4 
hours per configuration (the maximum possible taking into 
account the elevation restrictions for this source when ob- 
served from the northern hemisphere) and recovers emission 
from all spatial structures with angular scales between a few 
arcsec and 15 arcmin. 

All data were processed using the MIRIAD software 
package. To avoid the diffraction effects produced by point 
sources present in the field, for each of the brightest sources 
we imaged a small region around, and the clean components 
were Fourier transformed and subtracted from the visibili- 
ties. The residual visibilities, containing all the source struc- 
ture except for the offending point sources were then imaged 
and the point sources were added back into the SN 1006 
image in the image plane. To recover flux density contribu- 
tion from structures on angular scales larger than 15 arcmin 
(which is important in this case since SN 1006 is ~30 arcmin 
in diameter) we added single dish observations acquired in 
2002 with the Parkes 64 m radiotelescope placed in Aus- 
tralia. Also, since the primary beam of the VLA (the half- 
power beamwidth of a single VLA antenna) at A20 cm is 32 
arcmin, comparable to the size of the source, a correction 
was applied to the interferometric data taking into account 
the attenuation introduced near the primary beam edge. 

The final image has a synthesized beam of 7". 7 x 4". 8, 



cept for ij>o < 30° in the quasiparallcl model. In the case of quasi- 
parallel injection, the angle <po = 30° roughly separates cases 
with a centrally peaked morphology {(po < 30°) and a bilateral 
one ((j)o > 30°). 



^ The Very Large Array of the National Radio Astronomy Obser- 
vatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Figure 2. Azimuthal variation of the radio surface brightness ^(gsmaxi V^) for different aspect angles (po- Numerical results (solid lines) 
are normalized to ^(gsmaxi '""/S). Calculations are done for b = 0, (3 = 0, s = 2.2; the plots however are almost the same for different 
values of b, /3 and s. The models of injection are quasi-parallel (a), isotropic (b) and quasi-perpendicular (c). ^Smax = 0.98 (a) and 
5Smax = 0.97 (b, c). Experimental data for SN 1006 are shown for region I (in blue) and II (in cyan). They are measured at 1.5 GHz 
between 12 and 14 arcmin from the centre and are normalized to the maximum value of l-cr errors. Azimuthal profiles given by the 
approximate formula ((Sjl are shown by dashed lines. 



position angle 8°. 3, and an rms noise of 1 x 10"'* Jy/beam. 
When combined with Parkes single-dish data, the total re- 
covered flux is S=14.9 Jy, in excellent agreement withprevi- 
ous estimates from Green's (2006) catalogue of SNRq The 
new image is presented in Fig. [3] The accuracy of the field 
point sources serves for comparison to appreciate the quality 
of the new image that resolves the SNR radio features down 
to the same flue spatial resolutiorlf]. 



3.2 Results 

To determine the MF orientation we only considered the SE 
half of SN1006 (regions I and II on Fig. (J} because this part 
is quite spherical and therefore is more appropriate for com- 
parison with the numerical results obtained in Sect. [2] for a 
SNR in uniform ISM. From the radio map of SN 1006, we 
extracted the radial brightness profiles (along radii of the 
SNR projection separated on Aip — 12 degrees). The ex- 
perimental radial distributions are subject to pixel-to-pixel 
variations. In order to lower the possibility of error due to 
fluctuations in observational data, we calculated the aver- 
ages of brightness and 1-a errors within 12 to 14 arcmin 
from the centre of SN 1006 (where the maximum in radial 
distribution of the surface brightness is located). 

Experimental data are compared with the theoretical 
results on Fig. [2] The estimated aspect angle (j>o differs much 
for the polar-cap and the equatorial-rim models of SN 1006. 
From the numerical simulations, the best-fltted aspect angle 
is (po = 70° ± 4.2° for isotropic injection, 0o = 64° ± 2.8° for 
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Figure 3. Radio image of SN 1006, in linear scale (galactic co- 
ordinates). Axes outline regions I and II used in our analysis. 

quasi-perpendicular injection and (po = 11° ± 0.8° for quasi- 
parallel injection. Considering an isotropic i njecti on and 
equatorial-rim model for SN 1006, iRevnolcQ (Il996l ) found 
a similar aspect angle, (j)o ~ 60°. 



^ |http : //mrao ■ cam. ac.uk /snrs/ 

^ After submission of our paper, a new rad io image of SN 1006 
has been reported by ICassam- Chenai et al] 12008). This image, 
based on new observations, has an angular resolution of ^ 6" X 
9", comparable to that achieved in this work. It has a better 
noise level (20 /ijy/beam); however, it lacks the large spatial scale 
structures that were recovered in our image with the addition of 
single dish data. 



4 DISCUSSION 

4.1 Basic processes determining the azimuthal 
brightness profiles 

We found that the azimuthal profiles of the radio surface 
brightness of the Sedov SNR are almost independent of the 
shapes of the distributions of relativistic electrons and mag- 
netic field downstream of the shock. Close to the shock, the 
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azimuthal variations of the radio brightness of the SNR may 
therefore be approximately described as (Appendix A): 



Sei'f) CX <j(9o,eff (<y2, 0o)) Ab (So.cff ( V', 0o )) ''^"^^ ''^ 

where 

cosOo.cff ~ cos sin 00 

for isotropic or quasi-perpendicular injection and 



cos Oo < 



: sm Lp sm (po 



(3) 



(4) 



(5) 



for quasi-parallel injections. (Eq. ((3} is not valid for a 
centrally-brightened SNR.) 

In order to compare the approximate azimuthal profiles 
given by this formula with the numerical results, we plotted 
them on Fig. [5] with dashed lines. This figure shows that Q 
can be used as an approximation for the azimuthal distri- 
butions of the radio surface brightness in SNRs evolving in 
uniform medium and uniform magnetic field. 

An important consequence of Q is that there are two 
basic processes which determine the azimuthal distribution 
of the radio surface brightness in SNR evolving in uniform 
ISM and uniform ISMF: electron injection and compression 
(and/or amplification) of ISMF on the shock. 



4.2 Modified shock 

The estimations in Sect. 13.21 were obtained for a shock com- 
pression ratio (T = 4. Modified shocks may be responsible for 
larger compression, a = 7. Modeling the situation out of the 
Bohm limit, we may use such compression factor together 
with ((2]). Essentially, the estimated aspect angle does not 
change in this case. Namely, it becomes (po = 65° ± 4.2°, 
64° ± 2.8° and 11° ± 0.8° for isotropic, quasi-perpendicular 
and quasi-parallel injection respectively. Thus, even if the 
shock is modified but not in the Bohm limit, the estimation 
for the aspect angle obtained under the assumption of an 
unmodified shock, may well be valid. 

Our model is based on a classical MHD approach where 
the parallel and perpendicular components of the ambi- 
ent magnetic field are compressed in different ways, a s it 
is given by (E]> (e. g. iKorobeinikov fc Karlikovl (|l96Cll ) or 
iKorobeinikovl ( 199ll , Chapter 7.4) for the self-similar solu- 
tion of the problem of the strong point explosion in a con- 
stant magnetic field). This approach might not be valid in 
a limit in which the turbulent magnetic field dominates, as 
in case of the very efficient non-linear particle acceleration 
consistent with Bohm diffusion. The quasi-parallel theory 
assumes in this case that the turbulence is produced ahead 
of the shock, not downstream. The compression of the (al- 
ready turbulent) m agnetic field then does not depend on th e 
original obliquit y (Volk ot ajj l2003l : iBerezhko eral]|2002h . 
[Rakows kTet al.1 |2008) argue that shocks of different initial 
obliquity subject to magnetic field amplification become per- 
pendicular immediately upstream. We may model such situ- 
ation assuming crB(^o) is a constant instead of using ^ and 
B{f, Go) = B{f,-K/2). 

The estimations of the aspect angle in this case are 58° ± 
2.8° and 11°±0.9° for quasi-perpendicular and quasi-parallel 
injection respectively. Isotropic injection produces constant 
azimuthal profiles in this case, i.e. the modeled SNR looks 
like a rim for any aspect angle. 




Figure 4. Surface brightness distribution (in linear scale) of SNR 
with the same parameters as on Fig. |2] The surface brightness 
distribution is calculated assuming the quasi-parallel model of 
injection and the aspect angle </<o = 11°. The component of ISMF 
in the plane of the sky is parallel to the horizontal axis. The levels 
of brightness are spaced in the same way as on the Fig[l] 



The estimated aspect angles in the assumption of the 
Bohm limit are close to those obtained for an unmodified 
shock, except for the isotropic injection. Isotropic injection, 
together with the assumption that the downstream magnetic 
field is independent of obliquity is not able to reproduce bi- 
lateral morphology of SNR. There is no azimuthal variation 
of the radio surface brightness in this case, in agreement 
with dSl). 



4.3 Injection efficiency and obliquity 

It is worth emphasizing that our analysis may have some 
implications for the model of injection efficiency. 

Our argument against the quasi-parallel injection is the 
morphology of SN 1006 it should have in uniform ISM and 
uniform ISMF. Fig. |4] shows an image of SN 1006 in case 
of the quasi-parallel injection (? oc cos'^ ©s) and aspect an- 
gle (j)o = 11°. Since the ambient magnetic field should be 
almost aligned with the line of sight and injections prefer 
parallel shocks ('polar caps' directed toward and away from 
observer), the brightness distribution of SN 1006 should be 
centrally brightened (with one or two radio 'eyes' within 
thermal X-ray rim), contrary to what is observed. 

Our argument seems therefore to disagree with polar- 
cap morphology, and favour a NW-SE orientation of inter- 
stell ar magnetic field around this SNR. 

iRothenfiug et all l|2004l) suggested a criterion, which is 
considered as an argument on behalf of the polar-cap model 
for SN 1006. Namely, if SN 1006 has an equatorial rim the 
observer should see some emission between the bright limbs. 
Numerically, the value of the parameter Rtt/s defined as a 
ratio between total power coming from the interior and that 
from the limbs should be R„/'i > 0.5 if SN 1006 is a barrel. 
The most likely explanation for the smal ler value of /a 
is th at the visible limbs are polar caps (|Rothenflug et aU 
The value of t his ratio is Rn/s ~ 0.7 in radio 



( Rothenflug et al.]|2004l '). Therefore, based on the R-j^/s cri- 
terion alone, the radio data itself cannot give preference nei- 
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ther to equatorial-rim nor to pola r-cap model of SN 1006 . 
However, < 0.3 in X-rays (|Rotlienflug et all |2004 ). 

prefering therefore the latter model. 

To this end, our argument against polar caps are in con- 
tradiction with i?^/3-criterion applied to X-ray data. Note, 
that the fl^/3-criterion is obtained for cylindrical source of 
isotropic emission. It would be interesting to see how devi- 
ations from this assumption may affect the criterion. Our 
models assume a uniform ISM and uniform ISMF. Could it 
be possible to reproduce the bilateral morphology (i.e. to 
obtain (jjo > 30°) with quasiparallel injection if one consid- 
ers a gradient in the ISMF? Assuming the contrast of the 
ISMF between the NE and SE regions can be determined 
by the relationship Bse/Sne — (Sse/Sne)^^^ , a ratio of 4 
for (j)o = 45°, and 20 for 60°, could make the azimuthal pro- 
file of radio emission comparable to the observed one. Fur- 
ther investigation, including multi-dimensional modeling of 
SN 1006, as recent VHE 7-ray data could help us to under- 
stand the nature of morphology of this SNR. 



5 CONCLUSIONS 

Radio images of SNRs can be a useful probe to determine 
the orientation of ISMF around SNRs. We develop a method 
for determination of the aspect angle between ISMF and the 
line of sight, through the comparison between radio obser- 
vations and model calculations. It is based on a property of 
azimuthal variation of radio surface brightness. Namely, the 
azimuthal profile is almost insensitive to the downstream 
evolution of magnetic field and relativistic electrons. In con- 
trast, it is sensitive to the obliquity dependences of both 
the injection efficiency and the compression (amplification) 
factor of ISMF. The obliquity angle, the azimuthal angle 
and the aspect angle are geometrically related. The simple 
expression ([3} together with certain assumptions about the 
model of injection (where isotropic or quasi-perpendicular 
injection appears to be sugge sted by a number of arguments 
(Fulbrig ht fc Revnoldsl IToOOh ) may be used in order to fit 
the observed azimuthal profile and to approximately esti- 
mate the aspect angle. 

We applied our method to determine the aspect an- 
gle in SN 1006. A new radio image of SN 1006 obtained 
from archival VLA and Parkes data is reported. The im- 
age has very good angular resolution and has contributions 
from all spatial scales recovered. The use of the numeri- 
cal modelling together with the new observations reported 
here allows us to find the best-fitting value of the aspect 
angle in case of unmodified shock: (po = 70° ± 4.2° for 
the case of isotropic injection, tfjo — 64° ± 2.8° for quasi- 
perpendicular injection (equatorial-belt model of SNR in 
both cases) and (po = 11° ± 0.8° for quasi-parallel injection 
(polar-cap model). The angles estimated under assumption 
of modified shock are quite close to those given above. 

There are some limitations from our results on the 
model of injection. In quasi-parallel and quasi-perpendicular 
models, aspect angles are expected (and found) to be similar 
in both modified and unmodified shocks because the obliq- 
uity variation of injection dominates the obliquity variation 
of ISMF compression/amplification. Rather fiat observed az- 
imuthal radio profile strongly prefers <j!>o < 30° if injection is 
assumed quasi-parallel (Fig. [2^), i.e. for polar-cap model of 



SN 1006. In this case, SN 1006 should be centrally bright- 
ened, contrary to what is observed. Next, if one assumes 
the shock in SN1006 is so strongly modified that compres- 
sion of ISMF is independent of obliquity (Bohm limit), then 
the only possible injection model is quasi-perpendicular, be- 
cause isotropic injection and constant Ab results in constant 
azimuthal profiles of the radio surface brightness, again con- 
trary to what is observed. 

Rejection of the quasi-parallel injection model in 
SN 1006 means that the initial ISMF is directed from SE 
to NW and SN 1006 has a barrel-shaped, rather than polar- 
cap, morphology. 

With our results, we come to a puzzling issue which 
should be investigated in the future : lRothenflug et al.l (|2004h 
_R^/3-criterion applied to the X-ray emission seems to ex- 
clude the equatorial model of SN 1006 while the analysis 
of the azimuthal radio profiles seems to be against of the 
polar-cap scenario. 
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APPENDIX A: APPROXIMATE FORMULA 
FOR AZIMUTHAL VARIATION OF THE 
RADIO SURFACE BRIGHTNESS IN SEDOV 
SNR 

Here we derive an approximate formula for azimuthal varia- 
tion of the radio surface brightness. This formula allows one 
to avoid detailed numerical simulations and may be useful 
in situations where the approximate estimation of the as- 
pect angle is reasonable. In addition, the formula allows us 
to have deeper insight in the main factors determining the 
azimuthal variation of the radio surface brightness in SNRs. 

The downstream distributions of K and B in a Sedov 
SNR in uniform ISM and uniform ISMF are 



K (X <;{eo)K{f), B cx AB{eo)B{f, Go) 



(Al) 



If one neglects the small differences in downstream distri- 
butions of the p arallel and perpendicular components of B 
l|Revnolddll99^ . Fig. 1), then 



B(f, Go) 



(A2) 



The obliquity angle Oo is different for each radial sec- 
tor of 3-D object. It is determined, for any position within 
SNR, by the set {ip,f / Q,cf>o)- Integration along the line of 
sight gathers information from different radial sectors, with 
different obliquities. Let us determine the 'effective' obliq- 
uity angle by the relation 



6o,cff (V, <t>o) = 6o(<y3, 1, 4>o)- 



(A3) 



Actually, So, off for a given azimuth equals to the obliquity 
angle for a sector with the same azimuth lying in the plane 
of the sky (i.e. in the plane being perpendicular to the line 
of sight and containing the center of SNR). Qo varies around 
0o,cff during integration along the line of sight. The closer 
Q to the edge of SNR projection the smaller the range for 
variation of So and more accurate is our approximation. 
(Actually, we used g corresponding to maximum in radial 
brightness distribution which happens rather close to the 
shock.) 

Let us consider the azimuthal profile of the radio bright- 
ness Sq at a given radius g from the centre of the SNR pro- 
jection. 

With the use of Bo, off, the azimuthal variation of the 
radio brightness for fixed g may approximately be written 
from ((T)) as 



Sg (X <;(eo,eff)^B(eo,off)' 



V r'^ - 



(A4) 



The integral in HA4|) is the same for any azimuthal angle ip. 
The azimuthal variation of the radio brightness is therefore 
approximately determined by 



Se{ip) cx ?(eo,cff ((^, 0o)) AB(6o,off('y3,0o)) 



(s+l)/2 



(A5) 



The relation between the azimuthal angle tp, the obliquity 
angle Oo,cff and the aspect angle (j>o is as simple as 



cos Oo,cff = sin ip sin < 



(A7) 



for quasi-parallel injections (the relations are different be- 
cause we define the azimuth angle v^Smax = 7r/2 at the posi- 
tion of the maximum brightness). 

The accuracy of the approximation (|A5|I is shown on 
Fig. (2] The azimuthal profiles is sensitive to g in quasi- 
parallel case for aspect angles less than about 30°, i.e. 
for SNR with centrally-brightened radio morphology (? cx 
cos^ 6o and, for small aspect angles, Oo — > 7r/2 on the per- 
ifery of SNR and thus ? — > there). Therefore, the for- 
mula (|A5|I does not give correct profiles in the case of quasi- 
parallel injection for 0o < 30°, unless g 1. 



cos Oo.cft = cos ip sin < 



(A6) 



